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Abstract 



We compute the 0{a s ) QCD corrections to charged triple vector boson production 
at a hadron collider, i.e. the processes pp — > ZZW^ + X and pp — * W^W^W^ + X. 
Intermediate Higgs boson exchange effects, spin correlations from leptonic vector boson 
decays, and off-shell contributions are all taken into account. Results are implemented 
in a fully flexible Monte Carlo program that allows for an easy customization of kine- 
matical cuts and variation of the factorization and renormalization scales. We analyze 
the dependence of the differential cross sections under scale variations and present 
distributions where the QCD corrections strongly modify the leading-order results. 



1 Introduction 



With the advent of data from the CERN Large Hadron Collider (LHC), phenomeno- 
logical studies and interpretation of the data will require precise theoretical predictions for 
both signal and background processes. The calculation of higher-order terms in the QCD 
perturbation series thus becomes an even more important issue than at present. 

Triple vector boson production processes are of particular interest because they are sen- 
sitive to quartic electroweak couplings and they are a Standard Model background for many 
new-physics searches, characterized by several leptons in the final state. Recently, the QCD 
corrections for pp -> W+W~Z + X, ZZZ + X, W + W~W + + X and ZZW + + X have 
appeared in the literature [1-3]. With .fT-factors ranging from 1.5 to 2 at the LHC and 
a strong phase-space dependence, they show a behavior which is similar to that found in 
di-boson production in hadronic collisions, where QCD corrections have been known for a 
long time [4-6]. Thus, these next-to- leading order (NLO) calculations need to be taken into 
account for every phenomenological study involving triple vector boson production processes 
at the LHC. However, since vector bosons are identified via their leptonic decay products, 
the calculations should include the leptonic decays. Furthermore, intermediate Higgs con- 
tributions are not negligible since they can enhance the cross section significantly and lead 
to dramatic changes in the shapes of distributions for certain observables. 

In this paper we compute the next-to-leading order QCD corrections to the four processes 
pp -> ZZW- + X, ZZW+ + X, W+W-W+ + X and W~W + W~+X, with subsequent decay 
of the vector bosons into final-state leptons. All spin correlations involved in vector boson 
decays, the effects due to intermediate Higgs boson exchange and off-shell contributions have 
been correctly taken into account. Two of the four processes, namely pp — > ZZW~ + X and 
pp — > W~W + W~ + X have been computed for the first time here. The other two processes 
have first been presented in Ref . [3] , albeit without leptonic decays and without Higgs boson 
exchange contributions. 

The results of our calculations have been implemented in a fully flexible Monte Carlo 
program, producing total and differential cross sections at NLO as well as Les Houches event 
files at tree level. In this paper, we will always refer to a pp collider, having LHC in mind. 
However, in the program VBFNLO, which will be publicly available in the near future [7], 
protons can easily be replaced by anti-protons with a simple change in an input file. 

Our paper is organized as follows: in Sec. [2J we discuss the organization of the calculation, 
the different contributions to the leading order (LO) and NLO cross section and we describe 
the checks which we have performed. In Sec. [3J, results are presented for charged triple boson 
production at the LHC. We discuss the renormalization- and factorization-scale dependence 
and further show some sample distributions with strongly phase-space dependent i^-factors. 
Finally, in Sec. HI we give our conclusions. 
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2 Calculational details 



The calculation for the processes presented in this paper has been performed in complete 
analogy to the calculation for W + W~Z + X, with leptonic decays, described in Ref. [2]. 
In the following, the different contributions to LO and NLO cross sections are discussed 
in detail. Furthermore, the tests which we have performed to check the consistency of our 
results are described. 



2.1 Tree-level contributions 



We have evaluated the full set of Feynman diagrams for the different final states 
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using the helicity amplitude method described in Ref. [8]. All fermion mass effects have been 
neglected. The indices on the lepton pairs indicate that different generations are assumed 
for the decay products of the three vector bosons, i.e. interference terms due to identical 
leptons in the final state have been neglected. At LO, there are 209 diagrams for ZZW 
production and 85 diagrams for WWW production. 



a) 







Figure 1: A selection of Feynman diagrams for tree-level ZZW production. 



The tree-level diagrams can be grouped into three distinct topologies: 



a) The first one (see Fig. [Th) contains all diagrams where there is only one vector boson 
attached to the quark line, decaying further into 6 leptons. 



3 



b) The second topology (see Fig. [T]d) comprises all diagrams where exactly two vector 
bosons are attached to the quark line and then decay into two and four leptons, re- 
spectively. 

c) The third one (see Fig. [lb) consists of all diagrams where all the three vector bosons 
are attached to the quark line. 

These topologies give rise to different one-loop contributions as will be discussed later. 

The parts of the Feynman diagrams that describe the vector bosons decaying into leptons 
can be seen as effective polarization vectors, computed only once and used for different quark 
flavor flow. For example, for ZZW + production, the two different subprocesses ud — > 5£ + u 
and sc — > 5 i + v (with an anti-quark, s, having the same momentum as the up-quark in the 
first case, and identical d and c momenta) share the same effective polarization vectors. These 
effective polarization vectors are computed numerically at the beginning of the evaluation 
of the full matrix elements, at a given phase space point, and reused wherever they appear. 
This reduces the amount of time spent in the calculation of the tree-level matrix elements 
significantly. 

In the calculation of leptonic tensors, special care has to be taken in the treatment of 
finite-width effects in massive vector boson propagators. In our code, we use the modified 
complex-mass scheme as implemented in MadGraph, that is we globally replace m v with 
my — imyTv, while keeping a real value for sia 2 9w [14]. 

When intermediate Higgs boson exchange effects are included, particular care is needed 
in the generation of the phase space, since, for small Higgs boson masses (100-300 GeV), 
the Higgs boson width is very narrow. A Breit-Wigner mapping is needed for the efficient 
generation of this resonance. In the case of WWW production there are two distinct W + W~ 
pairs which can be produced from Higgs boson decay while the Higgs resonance appears only 
once in WW Z and ZZW production. For WWW production we, therefore, have generated 
the 1 — > 3 boson phase space using Dalitz plot variables which allow for simultaneous Breit- 
Wigner mappings of two different invariant di-boson masses and thus the two different Higgs 
resonances. This procedure is very important for good Monte Carlo statistics since the Higgs 
boson contributions can enhance the LO W + W~W + + X production cross section by up to 
a factor of 5 and the NLO W + W~W + + X production cross section by almost a factor of 4 
as shown in Table [TJ 

2.2 Real- emission contributions 

The total NLO cross section is given by the sum of the real emission contributions, the 
virtual contributions and a collinear term, a finite remnant after the initial-state collinear 
singularities are absorbed into the parton distribution functions (pdfs). The real emission 
and virtual contributions are separately infrared divergent in D = 4 dimensions. In order to 
deal with these divergences, dimensional reduction (D = 4 — 2e) has been used and we apply 
the dipole subtraction scheme in the formalism proposed by Catani and Seymour [9] . Exact 



4 



Higgs boson mass [GeV] 


a LO [fb] 


a NLO [fb] 


60 


0.1133 ±0.0002 


0.2141 ±0.0003 


120 


0.2256 ±0.0002 


0.3589 ± 0.0004 


160 


0.5964 ±0.0010 


0.8360 ±0.0016 


180 


0.4553 ± 0.0007 


0.6568 ± 0.0009 



Table 1: LO and NLO cross sections for the process pp — > l{ z/£ 2 z/£ 3 £3 +X at the LHC, 
within the cuts of Eq. ( 13. 3j) . for scales fj, F = \i R = 3m w and for four different Higgs boson 
masses. The quoted uncertainties on the cross sections represent Monte Carlo statistical 
errors only. 
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Figure 2: Sample of real-emission diagrams for ZZW production. 

expressions for the dipoles as well as for the finite collinear remnants have already been 
presented in Ref. [3]. The real-emission matrix elements can be divided into two different 
classes: 

a) Diagrams where the emitted particle is a final-state gluon (see Fig. (2^,, where the 
crosses represent possible gluon- vertex insertions) 

b) Diagrams where the emitted particle is a final-state quark, and a gluon is present in 
the initial state (see Fig. [2Jd). These diagrams can be obtained by a simple crossing of 
the diagrams of the previous class. As will be shown in Sec. EH these diagrams show a 
stronger scale dependence and are responsible for large .fT-factors in many distributions. 

The pre-calculation of the leptonic tensors, as effective polarization vectors, already applied 
for the LO matrix elements, leads here to an even larger increase in the evaluation speed. 
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2.3 Virtual contributions 




Figure 3: The three one-loop topologies appearing in the calculation of the virtual contribu- 
tions. 

One-loop corrections to the tree-level diagrams of Fig. [I] can be organized according to 
the three topologies encountered at tree-level: 



a) One-loop corrections to the diagrams with a single weak boson attached to the quark 
line, as in Fig. UK, give rise to simple vertex corrections, as illustrated in the corre- 
sponding Fig. [3h- This type of corrections exactly factorizes on the corresponding Born 
amplitude. 

b) Virtual corrections to the diagrams with topology as in Fig. [Tb give rise to the virtual 
diagrams illustrated in Fig. [3b- The sum of the four virtual contributions along a quark 
line will be called "boxline contribution" in the following. 

c) One-loop corrections to the diagrams with topology as in Fig. [lb give rise to the most 
complicated topology, Fig. [3b, where self-energy, vertex, box and pentagon corrections 
appear. The sum of all the virtual contributions along a quark line will be called 
"pentline contribution" in the following. 



Since there are only three colored partons in the real-emission diagrams, all infrared singular- 
ities appearing in the virtual contribution factorize on the Born amplitude. In conventional 
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dimensional regularization the virtual amplitude is then given by Q 



M v = M v + ^-C F [ ) T(l + e) 



,2 



2 3 „ 4tt 



e 2 e 



+ 



M B , (2.5) 



where Mb is the Born amplitude, s the square of the partonic center-of-mass energy and 
My is the finite contribution from the sum of all the one-loop diagrams. 

The boxline and pentline contributions have essentially the same analytic expressions 
found in the calculation of QCD NLO corrections in vector boson fusion processes, qq — > Vqq 
and qq — ► VVqq, discussed in Refs. [10] and [11] respectively, apart from crossing a final-state 
quark to the initial state, and performing then an analytic continuation. 

To deal with the finite boxline contribution, we have used the results obtained by a slightly 
modified version of the boxline routine discussed in Ref. [10]. This routine implements the 
Passarino-Veltman tensor reduction [12] and leads to quite stable results. 

The pentline reduction needs a more stable reduction procedure. We have implemented 
the method proposed by Denner and Dittmaier [13]. In addition, we have implemented a new 
calculation of the pentline contributions which reuses intermediate results for different vector 
boson polarizations. We have checked these results with the pentline routines computed in 
Ref. [11], after crossing and analytic continuation. The new pentline subroutines turn out to 
be 4.5 times faster and numerically more stable than the old code. For non-exceptional phase- 
space points, we found agreement for the two different codes at the 10~ 8 level. However, 
even with the increase in speed, this part of the code is still quite slow. Therefore, we have 
applied a trick, already used in Ref. [2] , to reduce the contribution of the pentagon diagrams: 
we have split the effective polarization vector e v of a vector boson of momentum qy into a 
term proportional to the momentum itself and a remainder e v 



e v = xyq v + e v . (2.6) 

The contraction of the pentline contribution with the component aligned along q v reduces the 
pentline itself to the difference of boxline contributions. Therefore it is possible to shift part 
of the pentline contribution to the less time-consuming boxline contributions and calculate 
the remaining smaller pentline contribution (the one obtained with the contraction with e v ) 
with less statistics, without changing the overall Monte Carlo statistical error of the total 
NLO result [11]. In practice we have chosen 

/ , x n e v ■ (qz 1 + qz 2 ) fnr? , 

tv ■ {qz 1 + qz 2 ) =0 => x v = ■ 2.7 

Qv • {Qz 1 + qz 2 ) 

for ZZW production and 

/ , , n n e v ■ (qw! + qw 2 + qw 3 ) , n c n 
ev ■ {qwj + qw 2 + qw 3 ) = o => x v = ■ ■ (2.8) 



qv ■ (qwi + qw 2 + qw 3 



for the WWW case. 



1 Using dimensional reduction instead, one needs to replace — 8Mb by —7 Mb, but this difference is exactly 
compensated by an analogous replacement in < 1(e) > , the integral of the real-emission counter-term. 
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2.4 Checks 



For all the triple boson production processes, we have performed numerous checks on the final 
results. All matrix elements in the LO and in the real-emission calculation have been checked 
individually against MadGraph and agree at the 10 -15 level. In addition, we have compared 
the LO cross sections against MadEvent [14] and HELAC [15] and find agreement within 
the statistical accuracy of the Monte Carlo runs (0.5% for HELAC and 1% for MadEvent). 
Furthermore, we have implemented Ward identity tests for the virtual contributions and 
checked the cancellation of divergences in the real emission against their counter-terms, as 
given by Catani and Seymour [9]. 

As a final and very important test, we have made a comparison with the already pub- 
lished results for the production of on-shell gauge bosons without leptonic decays of Ref . [3] . 
Since the authors of this paper have not included Higgs boson exchange and leptonic spin 
correlations in their calculation, we have neglected these contribution too, i.e. we have ne- 
glected the Feynman graphs with Higgs boson exchange and non-resonant contributions and 
we have used the narrow- width approximation for vector boson decay. In Tables [2] and 0, 



Scale 


program 


a LO [fb] 


a NLO [fb] 


1/2 x (3 m z ) 


VBFNLO 

Ref. [3] 


20.42 ±0.03 
20.2 ±0.1 


43.02 ±0.08 
43.0 ±0.2 


2m z + m w 


VBFNLO 
Ref. [3] 


20.30 ±0.03 
20.2 ±0.1 


39.87 ±0.08 
40.4 ±0.2 


(3 m z ) 


VBFNLO 

Ref. [3] 


20.24 ±0.03 
20.0 ±0.1 


39.86 ±0.07 
39.7 ±0.2 


2 x (3 mz) 


VBFNLO 
Ref. [3] 


20.03 ±0.03 
19.7 ±0.1 


37.39 ±0.07 
37.8 ±0.2 



Table 2: Comparison between our results and the ones of Ref. [3] for ZZW + +X production. 
All parameters and settings are taken from Ref. [3]. 

we show the comparison between the two sets of results, for different factorization- and 
renormalization-scale choices, here taken to be equal. Our NLO results agree at the 1% 
level, which is satisfactory, given the same level of agreement for the LO cross sections and 
the size of the Monte Carlo errors. 



3 Results 

The calculations described in the previous section have been implemented in a fully 
flexible parton-level Monte Carlo program, VBFNLO, which originally was developed for the 
prediction of NLO QCD corrections to vector boson fusion processes. The various triple 
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Scalp 


program 


a LO [fbl 


a NLO ffbl 


1/2 X (3 77lz) 


VBFNLO 

Rpf rn 


82.7 ±0.1 


152.5 ±0.3 

1 W o _|_ n c 




VBFNLO 

Ref. [3] 


82.8 ±0.1 
82.5 ±0.5 


145.2 ± 0.3 
146.2 ± 0.6 


(3 m z ) 


VBFNLO 

Ref. [3] 


82.8 ±0.1 
81.4 ±0.5 


143.8 ±0.3 
144.5 ± 0.6 


2 x (3 m z ) 


VBFNLO 

Ref. [3] 


82.4 ±0.1 
81.8 ±0.5 


136.8 ±0.3 
139.1 ±0.6 



Table 3: Comparison between our results and the ones of Ref. [3] for W + W W + + X 
production. All parameters and settings are taken from Ref. [3]. 

vector boson production options will be made publicly available soon [7]. The program 
allows for the calculation of cross sections and distributions in either pp, pp or pp collisions 
of arbitrary center of mass energy. In the following, we present results on ZZW and WWW 
production at the LHC, i.e. for pp collisions at ^fs = 14 TeV. 

The default electroweak parameters used in all plots are 

m w = 80.419 GeV m z = 91.188 GeV 

G F = 1.16639 x 10~ 5 GeV" 2 m H = 120 GeV. (3.1) 

Two other variables, a -1 = 132.507 and sin 2 (9 W ) = 0.22225, are calculated in the pro- 
gram using LO electroweak relations. We have used the CTEQ6M parton distribution with 
a s (mz) = 0.118 at NLO and CTEQ6L1 for the LO calculation [16]. All fermions are treated 
as massless and we do not consider contributions involving bottom and top quarks. The 
CKM matrix is approximated by a unit matrix throughout. The W- and Z-boson widths 
have been calculated in the program via tree-level formulas with one loop QCD corrections 
for the hadronic widths. For Higgs boson decays, approximate formulas are used which in- 
corporate running bottom-quark mass effects and off-shell effects in Higgs decays to weak 
bosons. For m H = 120 GeV, the resulting widths are 

T w = 2.0994 GeV , T z = 2.5096 GeV , T H = 0.004411 GeV . (3.2) 

In our calculations, the full leptonic final state is available and hence we determine cross 
sections for realistic acceptance cuts on the leptons. For WWW production, only a cut on the 
transverse momentum and the rapidity of the final-state charged leptons has been applied. 
For ZZW production we require, in addition, that the invariant mass of any combination 
of two charged leptons, m«, to be larger than 15 GeV, in order to avoid virtual-photon 
singularities in 7* — > £ + 1~ at low m u . Specifically, we require 

p Te > 10 GeV, \y e \ < 2.5, m u > 15 GeV (only for ZZW). (3.3) 
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solid: n F = n R = % 3 m w 
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Figure 4: Scale dependence of the LO and NLO cross section for 5 charged lepton final 
states within the cuts of Eq. A3.3\) . Left panel: variation of the renormalization and/or the 
factorization scale for ZZW + + X production. Right panel: same as in the left panel but 
for ZZW~ + X production. 



All results given below have been calculated for three different lepton families in the final 
state, i.e. interference terms due to identical particles have been neglected. Phenomenologi- 
cally more interesting are the cases of final states with electrons and/or muons. Considering 
decays of the three vector bosons into two generations of leptons each, the results for three 
distinct generations need to be multiplied by a combinatorial factor of four. This takes into 
account the presence of two identical vector bosons {ZZ and W ± W ± , respectively) and the 
corresponding symmetry factor of 1/2 which would appear when considering on-shell weak 
boson production. These factors are included in all figures. 

In Figs. H] and 0, the factorization- (/xp) and renormalization-scale dependence 
of the LO and the total NLO cross section is shown for all the different processes under 
investigation. At LO, there is no renormalization-scale dependence, since triple vector boson 
production is a purely electroweak process. Therefore, the scale variation is only due to the 
variation of the factorization scale in the parton distribution functions. The small variation 
at LO can thus be explained by the fact that the pdfs are determined in a Feynman-x 
range of small factorization scale dependence. At NLO, the dependence on the scales is 
more complicated. Since the factorization-scale dependence is quite small, the dependence 
under variation of // = //r = pip is almost completely dominated by the dependence on the 
renormalization scale and shows the expected a s (/x) dependence, i.e. the bigger the reference 
scale, the smaller the scale dependence. 

For a more detailed analysis, the different contributions to the total NLO cross section are 
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Figure 5: Scale dependence of the LO and NLO cross section for 3 charged lepton final 
states within the cuts of Eq. A3.3\) . Left panel: variation of the renormalization and/or the 
factorization scale for W + W~W + + X production. Right panel: same as in the left panel 
but for W~W + W~ + X production. 



shown in Fig. [6], for the example of ZZW + + X production. A qualitatively similar behavior 
is found for all triple vector boson processes investigated here. In the left panel of Fig. [6], 
the finite part of the virtual contributions (the My term in Eq. (I2.5p ). combined with the 
Born squared terms (including the LO contribution), show a remarkably small dependence 
under simultaneous variation of the renormalization and the factorization scale. This can 
be understood by a comparison with the factorization-scale induced LO variation given in 
Fig. |H Under variation of /ip, the virtual contribution shows the same behavior as the LO 
cross section, which means that the cross section decreases for small scales. Under variation 
of hr, on the other hand, the finite part of the virtual contribution increases for small scales, 
due to the increase in a s . These two opposing behaviors cancel to some extent and lead to 
the observed curve in Fig. [6j 

For the subtracted real-emission contributions and the finite collinear remnants, the ana- 
lysis of the scale dependence is somewhat more involved, since the finite collinear remnants 
depend non-trivially both on the factorization and on the renormalization scale. Moreover, 
these contributions include gluon-induced subprocesses like ug —> ZZW + d in addition to the 
quark-induced ones such as ud — > ZZW + g. In the right panel of Fig. [61 the real-emission 
contributions and the finite collinear remnants are therefore separately shown for each of 
these classes of subprocesses. 

In the finite collinear remnants, the quark- and gluon-induced contributions show op- 
posite behavior under variation of the scale. Due to these cancellations, the resulting scale 
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Figure 6: Scale dependence of the different contributions of the NLO cross section for pp 
ZZW + + X production at the LUC within the cuts of Eq. / EOl) . 



dependence and the size of their overall contribution is very small. The real-emission contri- 
butions arising from the quark-induced subprocesses show a similar scale dependence and are 
almost constant for the scales shown here. A comparatively large scale variation is observed 
in the real-emission terms of the gluon-induced contributions. These are also responsible 
for the large scale dependence of the overall real-emission term in the left panel of Fig. [6l 
This is not surprising since gluon-initiated subprocesses open up for the first time at NLO, 
and therefore, a LO-type scale dependence is expected. Gluon-induced subprocesses are also 
responsible for a large fraction of the i^-factor. For instance, the i^-factor for ZZW + produc- 
tion at fip = fin = 3mw is 2.1, whereas the i^-factor without gluon-initiated subprocesses 
only amounts to 1.5. 

In our analysis, we have also checked other scale choices, such as the invariant 3-vector 
boson mass or the minimal of the three vector bosons. We could not find an improved 
scale dependence, however, either in the cross section or in the distributions. This again 
can be understood since the dominating scale dependence comes from the gluon-induced 
subprocesses, which have to be considered as LO processes. 

For all processes studied, we have found a strong phase-space dependence of the size 
of NLO corrections. Thus, differential .fT-factors, defined as the ratio of NLO over LO 
differential distributions, 

can show a considerable variation. In the left panel of Fig. [TJ for instance, the invariant 
ZZ mass distribution in ZZW + + X production is shown for a Higgs boson mass of rrtH = 
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Figure 7: Differential cross section for the ZZ invariant mass in ZZW + + X production at 
the LHC. The Higgs boson mass used in the plot is m H = 170 GeV while u. F = u. R = 3 mw 
The ratio of the two distributions, defining the differential K-factor as given in Eq. Ii3.4\) . is 
shown in the right-hand panel. 



170 GeV. Here the Higgs boson contribution gives rise to the narrow peak at about mzz 
1 70 GeV. At tree level, the only Feynman graph with a Higgs boson exchange is the one 
depicted in Fig. [Th, where the Higgs boson decays into two Z bosons. This graph dominates 
near Mzz = n^H- In the right panel, we have plotted the differential i^-factor. Since the 
QCD corrections to the Higgs boson-contribution itself increase the cross section only by 
about 30% [17], there is a pronounced dip in the differential if-factor at about the Higgs 
boson mass. 

In Figs. [8] and [9], we give two more examples of the phase-space dependence of the NLO 
corrections for Vy + H /r_ Vy + + X production at m# = 120 GeV. In Fig. [HI the transverse- 
momentum distribution and the i^-factor of the highest-p^ charged lepton are shown while 
Fig. [9] shows the same for the charged lepton of lowest-pr- Variation of the if-factor up to 
70% is observed for the highest-pr lepton while for the lowest-px lepton we have variations 
up to 30% when considering "inclusive" event samples. This pr-dependence of the .fT-factors 
can be traced to the kinematics of the real-emission contributions. The rise is mostly due 
to events with high jets which are recoiling against the leptons. Imposing a veto on jets 
with px > 50 GeV leads to a fairly flat if-factor which, in addition, is close to unity for the 
lepton p T distributions (curves labeled "with jet veto" in Figs. [8] and |9]). Similar effects had 
previously been observed for vector boson pair production at the LHC [5]. 
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40 80 120 160 200 40 80 120 160 200 



Pt [GeV] p T [GeV] 

r 1 lep, max r 1 lep, max 

Figure 8: Differential cross section for the highest-p? lepton for \ir = pip = 3mw in 
W + W~W + + X production at the LHC. In the right-hand panel, the differential K-factors, 
as defined in Eq. ( 13.41) . are shown for inclusive events without jet cuts and also for a veto on 
jets with pr, jet > 50 GeV. 

4 Conclusions 

The simulation of triple vector boson production at the LHC is important for two reasons. 
These processes are a Standard Model background for new-physics searches which are char- 
acterized by multi-lepton final states, and secondly they are sensitive to quartic electroweak 
couplings. In this paper, we have presented first results for the full NLO differential cross 
sections for WWW and ZZW production, with all spin correlations from leptonic vector 
boson decays, intermediate Higgs boson-exchange effects and off-shell contributions taken 
into account. Results are collected in a fully flexible Monte Carlo program, VBFNL0 [7]. 

When varying the factorization and the renormalization scale p = pp = p^ up and down 
by a factor of 2 around the reference scale p = 3mw, we have found a scale dependence 
of about 5% for the LO cross section and of somewhat less than 10% for the NLO cross 
section, for WWW production. For the ZZW case, the LO scale dependence is around 
1%, whereas the dependence of the NLO cross section is around 13%. These variations are 
in the expected range for the NLO scale dependence, while the LO variations have to be 
considered anomalously small, due to the absence of initial-state gluon-induced subprocesses. 
The large A'-factors (of order 2 and even larger in some phase-space regions) demonstrate 
the importance of including the NLO QCD corrections on top of the LO predictions. 

The differential i^-factors for several distributions for both of these processes are highly 
dependent on the Higgs boson mass. In general we observe that the larger the contributions 
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Figure 9: Differential cross section for the lowest-px lepton for p,p = u.r = ?>my/ in 
W + W~W + + X production at the LHC. In the right-hand panel, the differential K-factors, 
as defined in Eq. Ii3.4\) . are shown for inclusive events without jet cuts and also for a veto on 
jets with Pt, jet > 50 GeV. 



from the Higgs boson are, the smaller is the i^-factor. In the case of the W + W~W + + X 
production, with p,p = hr = 3mw, for example, the K-factor decreases from 1.6 for a Higgs 
boson mass of 120 GeV to 1.4 for a Higgs boson mass of 150 GeV. At the same time the LO 
cross section increases by more than a factor of 2. Therefore, in all simulations, the Higgs 
boson contribution has to be taken into account in order to obtain a valid prediction for 
the cross sections and the ^-factors. Besides these large .fT-factors, we have also found a 
strong phase-space dependence of the size of the NLO corrections which shows that a mere 
multiplication of distributions by an overall .ff -factor is not sufficient. 
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